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Abstract 

We propose and develop a formalism to describe and constrain statistically anisotropic primor- 
dial perturbations. Starting from a decomposition of the primordial power spectrum in spherical 
harmonics, we find how the temperature fluctuations observed in the CMB sky are directly related 
to the coefficients in this harmonic expansion. Although the angular power spectrum does not 
discriminate between statistically isotropic and anisotropic perturbations, it is possible to define 
analogous quadratic estimators that are direct measures of statistical anisotropy. As a simple 
illustration of our formalism we test for the existence of a preferred direction in the primordial per- 
turbations using full-sky CMB maps. We do not find significant evidence supporting the existence 
of a dipole component in the primordial spectrum. 
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I. INTRODUCTION 



At the time Einstein decided to apply his equations to the universe, observational data in 
cosmology was rather scarce, if not inexistent. Hence, instead of relying on observations to 
constrain the spacetime metric, he postulated the "cosmological principle," the isotropy and 
homogeneity of the universe. Similarly, guided by theoretical prejudice, Einstein assumed 
that the universe was static, which, incidentally, is what forced him to introduce a cosmo- 
logical constant in his equations. It turns out that the universe is actually homogeneous 
and isotropic on large scales, as confirmed, for example, by the homogeneity of the distribu- 
tion of luminous red galaxies Q and the isotropy of the cosmic microwave background Q] • 
However, the universe is not static, and Einstein's staticity assumption prevented him from 
predicting the expansion of the universe later discovered by Hubble. 

The cosmological principle has a similar, but formally independent counterpart for cos- 
mological perturbations. It states that perturbations are statistically homogeneous and 
isotropic, that is, that their correlation functions are invariant under translations and rota- 
tions. These assumptions are so ingrained and integrated into our treatment of cosmological 
perturbations, that often we are unaware of them. As a consequence, the statistical ho- 
mogeneity and isotropy of cosmological perturbations has received little attention in the 
literature, and only recently the large scale anomalies observed in the cosmic microwave 
background radiation (CMB) have led to a spurt of interest in statistical anisotropy. 

In fact, the data provided by the WMAP cosmic microwave background experiment Q], 
shows hints of violations of statistical isotropy. These hints include the alignment of the 
quadrupole and octopole unlikely correlations between "multipole vectors" |4], a north- 
south asymmetry in the CMB sky [5( and evidence for the existence of a symmetry plane 

E|. On the other hand, a measure of statistical anisotropy proposed in j^] shows no 
evidence for statistical anisotropy, and it has been also suggested that the previous hints are 
artifacts of the heavily processed maps used to analyze the data p. 

In any case, regardless of whether the CMB sky is truly statistically anisotropic or not, the 
statistical isotropy of cosmological perturbations should be subject to empirical test, rather 
than taken on faith. Whereas most studies have focused on the statistical isotropy of the 
CMB sky, in this work we analyze the isotropy of the primordial perturbations themselves. 
Verifying whether cosmological perturbations are statistically anisotropic is a way to confirm 
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whether our basic assumptions about structure formation, and ultimately our universe, are 
correct. It provides a novel way to test our current understanding of the origin of structure, 
embodied in inflationary models, and it might provide information about the universe at 
high energies and long wavelengths. 

In this work we are not concerned with the mechanism responsible for the generation of 
statistically anisotropic primordial perturbations, although if it were impossible to seed them 
our analysis would be perhaps an academic exercise. The most straightforward mechanism 
to explain statistical anisotropies involves isotropy violations in the background spacetime 
itself. These departures from isotropy can arise from non-trivial spatial topologies Q| or 
departures from the background Friedmann- Robert son- Walker metric However, these 

explanations require the breaking of isotropy from the onset, and they also clash with our 
understanding of the origin of structure. Alternatively, statistical anisotropies might arise 
from coherent magnetic fields in the universe [12|. If these magnetic fields can be regarded 
as perturbations in an isotropic spacetime, their contribution to the metric perturbations is 
quadratic in the fields, and hence presumably highly suppressed. On the other hand, if they 
significantly contribute to the energy density of the universe, one expects deviations from 
isotropy at the level of the background spacetime itself, violating again isotropy from the 
onset. In this paper we assume that the universe is homogeneous, isotropic, spatially flat and 
practically infinite. But even homogeneous and isotropic universes can be dominated by sets 
of non- vanishing non-scalar background fields if the latter appear in special configurations, 
like for example, vector field triads. In these universes, perturbations are expected to be 
statistically anisotropic [l3j . In general the existence of these non- vanishing fields signals the 
spontaneous breaking of Lorentz-invariance ^| (rotational invariance to be more precise), 
so the study of statistical anisotropies in the primordial perturbations opens a new window 
to test the basic symmetries of nature. 

The paper is organized as follows. In Section |H] we precisely define what we mean by 
statistically isotropic or anisotropic primordial perturbations. It is in this section where we 
introduce a quantitative characterization of primordial statistical anisotropies. In Section UTTl 
we connect statistical anisotropies in the primordial perturbations to statistical anisotropies 
in the cosmic microwave background temperature fluctuations. We define statistics that 
capture the amount of statistical anisotropy in the CMB and at the same time mirror the 
amount of anisotropy in the primordial density ripples, and we also discuss the relation 



3 



between our statistics and the bipolar spectrum of Hajian and Souradeep |8|. Section llVl 
contains a description of strategies to test and constrain the statistical isotropy of the pri- 
mordial perturbations. This task is simplified if there is a cut-off in the anisotropies at a 
particular arbitrary multipole, which is the case we mainly consider. As an illustration of the 
formalism, we look for the simplest departure from statistical isotropy: a preferred direction 
in the primordial perturbations. We apply our statistics to the full sky CMB maps and do 
not find evidence of statistical anisotropy. Finally, in Section IS we draw our conclusions. 
We have also included an Appendix where we summarize most of the formulae that we need 
in our derivations. 



II. STATISTICALLY ANISOTROPIC RANDOM FIELDS 
A. Random Fields on Euclidean space 

Consider a real Gaussian random field in Euclidean (flat) three-dimensional space, 
At this point it is not important to know what $ actu ally is, although the reader can think 
of $ as the Newtonian potential (in longitudinal gauge 
space consist of translations, rotations and reflections 
how random fields behave under these transformations. The Gaussian random field <3> is 
statistically homogeneous if its momenta are invariant under translations, 



LllV is, < 

loi llq ) . The isometries of Euclidean 
17| . It is therefore useful to analyze 



($(£)) = ($(£+*)), ($*(£)$(#)) = ($*(£ + t)$(y + t)) VfGM 3 , (1) 

where ( ) denotes ensemble average. 1 It is not essential to assume that the field is Gaussian. 
This assumption just allows us to concentrate on the average (&(x)) and the correlation 
($*(x) $(?/)), which uniquely characterize the statistical properties of the random field in 
this case. A non-Gaussian field that does not satisfy equations (0) is not statistically homo- 
geneous either. 

The average ($(af)) of a statistically homogeneous random field is a constant in space, and 
hence is trivially invariant under rotations. A Gaussian random variable $ is statistically 



1 Because $ is real by assumption, $ = $*. We include the complex conjugation in the two-point function 
for convenience. 
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isotropic (around O) if its momenta are invariant under rotations (around O), 

($(£)) = md + K-(x-0)) (2) 
= + (f -0))$(0 + ft- (y-O))) VReSO{3). 

If we assume that the field is statistically homogeneous, equations (j2J) are satisfied for all 
points O G M 3 if they are satisfied for a single O (say, the origin O = 0.) Conversely, one can 
also show that if a random field $ fulfills equations (j2J) for all possible points O, then it also 
fulfills equations In other words, if a random field is statistically isotropic around all 
points, it is statistically homogeneous. Note that it is not necessary to discuss the statistical 
properties of <3> under inversions. The inversion of a single point can be accomplished by a 
translation, and the inversion of two points can be accomplished by a rotation. 

Often it is convenient to decompose spatial functions in eigenvectors of the translation 
operator. In flat space these are the plane waves exp(ik -x). Let us hence consider a Fourier 
decomposition of the random field, 

The components $(/c) are random fields in Fourier space. In this space, the requirement of 
statistical homogeneity translates into 

<$(£)) = %5 {3) (k) and ($*(£)$(#)) = (2n) 3 5 {3) (k- k')^-V(k), (4) 

where $o is the (constant) mean of $ and V(k) denotes the power spectrum, 2 which is an 
arbitrary positive definite function 3 of k. On the other hand, the requirement of statistical 
isotropy (j2J implies 

V(Rk) = V(k) VR. g 50(3) (5) 

that is, that the power spectrum only depends on the magnitude k of the wave vector k. 
Because any positive function can be written as a square, we can express the power spectrum 

as 

V{k) = k 3 $*{k)${k), (6) 



2 Note that our power spectrum has a slightly unconventional normalization. 

3 One can show this by considering the two-point function of a smoothed field & s : < ($>*(x)& s (x)}. 
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where $ is a real function by definition. Then, the random variable $ is statistically isotropic 
if and only if $ only depends on k. 

In this paper, we consider statistically anisotropic primordial perturbations, that is, pri- 
mordial perturbations whose power spectrum does not only depend on the magnitude k. It 
is hence going to be useful to decompose $(fc) into components that transform differently 
under rotations of the wave vector, 

&(k) = y/4ft^2$i m (k)Yi m (k), (7) 

Im 

where k = k/k and the Y\ m are spherical harmonics. Because $(fc) is real, the complex 
functions satisfy $* m = (— l) m $/_ m . This decomposition of $ carries over to the power 
spectrum itself. Inserting equation ((7J) into equation © and using relation (|A8|) we find 
that 

V(k) = ^^Vi m {k)Y lm {k), where (8) 

im 

V lm (k) = Yl (-ir^'^W^^Wi-mi^a.mal^m), (9) 

limi,! 2 m 2 

and the real coefficient D is given by equation (jA6|) . If the random field is statistically 
isotropic, V\ m = for I > 1. On the contrary, if $ is statistically anisotropic, there exists 
at least one I > 1 and m e {—/,...,/} such that V\ m ^ 0. For simplicity we shall express 
our results mostly in terms of Vi m , though the reader should be aware that the actual 
free parameters that characterize the statistical anisotropies are the $; m (fc). Observe that 
because V is real, V* m = (—l) m Vi- m , which also follows from the relation 0. 

B. Random Fields on a Sphere 

Sometimes, one is interested in random fields A(n) that are defined not on M 3 , but rather 
on a two-sphere, n 6 § 2 . The standard case is the temperature fluctuations on the cosmic 
microwave background sky. Because the isometries of a sphere consist only of rotations 
and inversions, for Gaussian fields it suffices to consider the properties of the random field 

n 

under rotations (the behavior of the CMB under inversions has been studied in A 
Gaussian random field is statistically isotropic if the average (A(n)) and the two-point 
function (A*(n)A(m)) are invariant under rotations, 

(A(n)> = (A(R,-n)), (A*(n)A(m)} = (A*(TZ ■ n)A(R, ■ m)) \I1Z G SO(3). (10) 
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In the applications we have in mind, random fields on a sphere receive weighted contribu- 
tions from an infinite number of Fourier modes. For instance, the t emp erature fluctuations 
measured by an observer at xq in the direction n can be written as jl6( 

where $ is a random field in Fourier space and A{k } n) is a real function. We shall specify 
what $ and A(k,n) represent below. At this point, we just want to know what constraints 
statistically isotropy of A T (x ,h) imposes on the quantities A(k,h). First, note that it 
follows from equation ([TT]) that if $ is statistically homogeneous, (At(xo>^)) is invariant 
under rotations if A(0,7?ii) = A(0,n). Furthermore, inserting the expansion (fTTj) into the 
second equation in (fTT)|) we arrive at the condition 

d 3 k 



J 



A(k,m)\ V(k) - \A{lZk,lZn)\ V(lZk) =0. (12) 



(2tt) 3 

Therefore, A^(o;o,n) is statistically isotropic if 

A{Kk,Kn) = A(k,n) and V(TLk) = V(k) VIZ G SO(3). (13) 

Note that the two-point function (A^(xo, n) A^(xo, rh)) does not depend on xq if the random 
field $ is statistically homogeneous. 



III. TEMPERATURE ANISOTROPIES 



Our main concern here are the contribution of scalar perturbations to the temperature 
fluctuations in the cosmic microwave background. Let us denote by At(x, h) the tempera- 
ture fluctuations measured by an observer at x looking at direction n in the sky. As hinted 
above, it is more convenient to study the Fourier transform of the previous fluctuations, 
Aj<(k, n). In linear perturbation theory, the evolution of At(/c,tz) is described by a set of 
decoupled differential equations for each mode k. The linearity of these equations allows 
us to separate the initial conditions from the evolution. In particular, for adiabatic initial 
conditions one can write 

A T (k,h) = A(jfe,n)$(jfe), (14) 

where A(k, h) is the temperature fluctuation one obtains by evolving the perturbation equa- 
tions from initial conditions where $(&;) = 1, and $(fc) is the value of the Newtonian po- 
tential at a sufficiently early time, well into the radiation dominated era. Note that this 
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split between evolution and initial conditions is arbitrary to some extent: what we consider 
as initial conditions here might be regarded for instance as the outcome of an inflationary 
stage. 

Recall now our discussion of statistical isotropy of a random field defined on a sphere. 
It follows from equation (|13|) that violations of statistical isotropy in the CMB sky might 
originate from two sources: from the initial conditions (encoded in the power spectrum 
V), or from the evolution (encoded in the transfer function A(k,h)). In this work we 
assume that the evolution does not generate statistically anisotropic perturbations, that is, 
we assume that A(k,n) = A(TZk,TZn) for all rotation matrices 1Z in SO(3). This is the 
case in cosmological models where the universe does not contain any non-scalar background 
field, like in the standard ACDM scenario. 4 As a consequence, the differential equations 
that describe the evolution of the perturbations are invariant under rotations of k and h, 
and it follows that Aj-(k, h) is a function of the two scalars one can construct out of k and 
n, namely k and k ■ n, A(k, h) = A(k, k ■ h). This property then allows us to expand A in 
Legendre polynomials Pi, 

A(k, ft) = ^(2Z + l)(-i) l Mk)Pi(k ■ n). (15) 
i 

This expansion will turn to be useful to compute the two-point function of the temperature 
anisotropies. 

A. Angular Two-point Function 

It is customary to characterize the properties of the temperature fluctuations by the 
angular two-point function (af m air m r), where the a\ m are the coefficients in the spherical 
harmonic decomposition 

A T (x ,n) = ^2ai m Yim(n), (16) 

lm 

and xq is our position in space (because of statistical homogeneity, the statistical properties 
of the temperature fluctuations do not depend on i .) We would like to express this angular 

4 In cosmological scenarios that contain non-scalar background fields, it is possible to construct scalar 
perturbations that are not invariant under rotations of k For example, if there exists a background 
vector field A in the universe, the scalar k ■ A $ is not invariant under k — > 1Z ■ k. 
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two-point function in terms of the initial power spectrum. Using the expansion (JT5j) and the 
addition theorem (jA2|) we find first 

/d 3 k 
j^y 3 Mk)<S>(k)YL(k). (17) 

Then, using the definition of the power spectrum (J1J and its expansion in spherical harmonics 
(JEJ) we arrive at our sought result 

( a z*imi a /2m 2 ) = {-i) h ~ h ^ D{k, mi] I, m\l 2 , m 2 )K(h, l 2 ; I, m), where (18) 

lm 

/dk 
—Al(k)A h (k)V lm (k), (19) 

and we have used equation (jA8|) in the Appendix. The structure of the two-point function 
can be understood in the language of addition of angular momenta. The l,m multipole of 
the primordial anisotropies can be thought of as a state with angular momentum quantum 
numbers I and m. The temperature perturbation a; imi "creates" angular momentum with 
quantum numbers h,mi_, and, similarly, a^ 2m2 destroys angular momentum with quantum 
numbers l 2 , m 2 . Therefore, the contribution of Vi m to the two-point function vanishes if the 
angular momenta l±, ni\ and I, m cannot add to l 2 , m 2 . In fact, D is zero unless m 2 = rrii + m 
and l 2 E {\h - /|, . . . ,/i + I}. 



B. Large Scale Anisotropies 

In general it is only possible to compute the transfer functions A; in equation ()19|) nu- 
merically. In order to provide a better handle on statistical anisotropies, in the following we 
consider a limit where it is possible to analyze them analytically, namely, on large angular 
scales. As it is well known, on these scales, neglecting the integrated Sachs- Wolfe effect, the 
Sachs- Wolfe relation 

A,(*0 = (20) 

hoMs Q. w„ re n, ienotes a. _ g dJL t „ « It „ ,„„ t0 

realize that this result is an approximation with a limited range of validity. In particular, 

5 By "last scattering" we mean the decoupling of photons and matter that occurred around redshift 
z = 1100. 
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at I < 8 the integrated Sachs- Wolfe effect modifies 6 the relation (|2Uj) an d at I > 30, the 
baryon velocity significantly contributes to the temperature anisotropics 19]. Therefore, the 
approximation (|20j) is only valid around 8 < / < 30. 

Just as one parametrizes the primordial spectrum in terms of an amplitude and a spectral 
index, let us parameterize the different components of the anisotropies $/ m as power laws. 
It follows then from equation (^J) that Vi m is a superposition of different power-laws, with 
different amplitudes and spectral indices. Let us assume for simplicity, and without loss of 
generality that 

k x 



Vim = Am ( -j^J ■ (21) 

In this case, the integral in equation (fTTIj) can be explicitly evaluated. Using equations 
(11.4.33) and (15.1.20) in Q we arrive at 



K(h, l 2 ; lm) 



2 n - 4 7rr[3-n]r[(Z 2 + /i + n - l)/2] A 



l in 



9 T [(h - l 2 + 4 - n)/2] T [(l 2 -h + A- n)/2) T [(Z 2 + h + 5 - n)/2)) (KD s ) n - 1 

(22) 



Hence, given the amplitude A and the spectral index in any given anisotropy multipole, 
one can exactly predict what the two-point function on large angular scales should be. In 
general however, we do not know what A or n, are. Under these circumstances, equation 
(I22|) in combination with (j!8)l could be used to fit for the values of A and n (see below.) 



C. Angular Power-Spectrum 

If primordial perturbations are statistically isotropic (Vi m = for / > 1), one can easily 
verify that the two point function ()18|) is diagonal. In this case it is customary to define the 
angular power spectrum Ci through 

( a l 1 m 1 a hm 2 ) = filil2$mim2- (23) 

But if perturbations are statistically anisotropic the two-point function is not diagonal, and 
the definition of the angular power spectrum above is meaningless. In any case, one never 
measures the CVs directly. Instead, one uses an appropriate estimator to infer their values. 

6 We assume a non-zero cosmological constant with Q\ sa 0.7. Since we do not know what is causing cosmic 
acceleration, this is just a guess. Other forms of dark energy lead to different ISW effects. 
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Here we shall define the CVs as the expectation value of the estimator that is commonly 
used to compute the latter, 

1 1 

Ci = (Q) = — J2 (*lflhn)- (24) 

m=—l 

The angular power spectrum defined this way is invariant under rotations. It reduces to 
()23|) if perturbations are statistically isotropic and it is also closer to the quantities that 
we actually measure. Note however that in the presence of statistical anisotropics, the 
CVs defined in this way do not completely characterize the statistical properties of the 
temperature fluctuations, since the two-point function is not diagonal. 

We want to find out now how the angular power spectrum is related to the primordial 
power spectrum. To this end, we substitute equation (fTSj) into the definition (J2U). Using 
then the identity equation (jA9|) we get 

/dk 
-\^(k)\ 2 V m {k). (25) 

Therefore, Poo ac ts like an effective power spectrum. If primordial perturbations are sta- 
tistically isotropic, then the effective power spectrum reduces to the conventional power 
spectrum, Poo(^) = "P{k). Otherwise, the effective power spectrum receives contributions 
from all anisotropy components, 

Poo = ]TA; 3 i>L<i /m , (26) 

l,m 

where we have used equations (|A3J) and (|A4|) . Let us stress that the angular power spectrum 
defined in equation (|24j) does not differentiate between the different anisotropy multipoles. A 
scale invariant monopole perturbation ($oo) leads to the same set of CVs as a scale invariant 



anisotropy in any other multipo 
CMBFAST |2l|, or CMBEASY 



e (i>; m ). Hence, conventional numerical codes, such as 
22) 1 can be used to determine the imprints of statistical 
anisotropies on the angular power spectrum. Alternatively, one can use equation (|2"2"|) to 
explicitly evaluate the C;'s on large — but not too large — angular scales. 

IV. LOOKING FOR PRIMORDIAL STATISTICAL ANISOTROPIES 

The power in the different multipoles of the primordial perturbations Vi m (k) uniquely 
determines the different two-point functions of the temperature anisotropies. The latter are 
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observables, provided that sufficient realizations of the random field are available. There- 
fore, one can hope to determine what the Vim are, or at least, place upper limits on their 
amplitude. This section is devoted to the ultimate goal of finding ways to measure V[ m . 

A. Determining K 

The first obstacle we have to find in our quest to measure V\ m is that the power is 
hidden behind a convolution, equation (fTTIJ) . Even if primordial perturbations are statisti cally 
isotropic it is difficult to invert that integral in order to find the primordial power [23]. 
Therefore, as a first step towards our goal we shall concentrate in determining the integral 
K in equation ()19|). Because K(li, I2; l,m) vanishes when Vi m = 0, we can think of if as a 
measure of the anisotropic power; a non-vanishing K for I 7^ would automatically imply 
the existence of statistically anisotropic primordial perturbations. The advantage of focusing 
on K is that we do not have to know what the transfer functions Ai(k) are. Therefore we 
can analyze all angular scales without having to make any assumptions about the underlying 
cosmology. 

The task of determining some of the values of K is greatly simplified when the anisotropies 
have a cut-off at at multipole I = L, $; m = for I > L. In this case, it follows from equation 
(JHJ) that the power V\ m is cut-off for multipoles I > 2L. Let us drop the factor of two and 
assume that Vi m = for I > L, with L even. In this case, K(l±, I2] I, m) vanishes for I > L 
and we can use the two-point function to directly determine the values of K for I = L. At 
the end of this subsection, we shall comment on how to iteratively establish whether there 
exists a putative cut-off at / = L. 

Consider equation (fTHj) with l\ = I and l 2 — l + L. Then, due to the properties of angular 
momentum addition 

(a* m a l+L<m+M ) = {-i) L D{l, m; L, M\l + L,m + M)K(l, l + L;L,M), (27) 

which allows to directly solve for K in terms of the two-point function. However, the two- 
point function is not directly observable. In order to measure K we need a sufficient number 
of realizations of the random process. At this point, note that the value of K does not 

7 Since we have access to only one universe, this condition is not automatically guaranteed. See the discus- 
sion of cosmic variance below. 
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depend on m. Hence we can define the following unbiased estimator, 

K(l, l + L;L,M) = 4— Y —, aLai + L > m+M -. (28) 

V ' ' ' ; 2/ + 1 „ D{l,m;L,M\l + L,m + M) V ; 

For L = M = 0, our estimator K gives the Cj's, equation (|24|) . Therefore, K can be regarded 
as a generalization of the angular power spectrum. However, note that for L ^ the variable 
K is not a scalar, and, in particular, does not transform under an irreducible representation 
of the rotation group. Let us stress here an additional property of K. Namely, because 
K is quadratic in the temperature fluctuations, our study of isotropy is decoupled of issues 
about non-Gaussianity. If the expectation value of K does not vanish for L > 0, primordial 
perturbations are statistically anisotropic, irrespective of whether temperature anisotropies 
are Gaussian or not. However, we cannot measure (K) directly, and the different realizations 
of the random process K are expected to fluctuate around its expectation value. It is here 
where we have to make additional assumptions about the nature of the fluctuations. In order 
to estimate the variance of K, let us assume for simplicity that perturbations are Gaussian 
and statistically isotropic. This assumption is justified because we are attempting to detect 
the existence of statistical anisotropies. Then, under these conditions, we find that the real 
and imaginary parts of K are uncorrelated, and for L ^ 

((K-m, , + L;L, M)) = ggg. g _ l _L _ - _ . (29) 

We encounter here the same cosmic variance problem one faces when trying to measure the 
angular power spectrum. Because for a given I we only have a finite number of samples of 
the two-point function, 21 + 1, we cannot measure (K) with infinite precision. By the way, 
note that if the value of the estimator happens to be consistent with zero, then we guessed 
incorrectly, and primordial anisotropies are not cut-off at I — L, but eventually at I = L — 2. 

B. A Different Estimator for K 

As we shall see, the orientation-dependence of our statistic K makes the analysis of 
temperature anisotropies rather cumbersome. Therefore, it may be advantageous to use an 
estimator of K that possesses well-defined properties under rotations. In order to construct 
it, let us consider the quadratic combination 

B(l 1} l 2 ;l,m) = a* imi a; 2m2 (-l) Wl (Z 1 ,-mi;i 2 m 2 |i,m), (30) 

mi,ra 2 
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where {Ii,mi,l2,m2\l,m) is a Clebsch-Gordan coefficient. The expectation ^value of B is 
nothin g e 
shown 



se but the bipolar spherical harmonic of Hajian and Souradeep [8]. It can be 
24j that under rotations B transforms under an irreducible representation of the 
rotation group, 

B(h, l 2 ; Z, m) - £ Dl mrn B{h, k, I, m), (31) 
in 

where D is the Wigner matrix of the rotation. Using equation ([180. we can express the 
expectation value of B (the bipolar coefficients) in terms of our quantity K, 

< B > = Yl ^(-0' 2 ^ 1 ^(^i,^i;^M|/ 2 ,m 2 )ir(/ 1 ,/ 2 ;L,M)(-ir i (/ 1 ,-m 1 ;/ 2 ,m 2 |/,m). 

(32) 

Again, the previous expression simplifies if there is an anisotropy cut-off at multipoles I > L. 
In that case, one can use the previous expression as a guidance to derive an unbiased 
estimator of K, 

i L B(l,l + L;L,M) 



K(l,l + L\L,M) 



£ m (-l) m £>(Z, m; L, M\l + L,m + M)(l, -m;l + L,m + M\L, M) ' 

(33) 



which is guaranteed to transform appropriately under rotations. However, as opposed to 
the previous case, we only have one realization of the quantity we are trying to estimate, so 
our estimator suffers from a larger cosmic variance. As a matter of fact, the mean square 
fluctuations of K are given by 

(£*)(., l + L;L,M))= m; L> m + , (34) 



where we have used equations (JA6[) and (jA7[) and, once again, we have assumed that pertur- 
bations are statistically anisotropic. It is apparent from the last equation, that the variance 
of K, equation ()29jl . is smaller that the one of /C. In this paper we shall concentrate on the 
lower variance estimator K; we reserve the study of /C and related quantities to a hopefully 
forthcoming publication. 

To conclude this section, let us comment on the relation between our integral K and 
the bipolar spectrum we briefly mentioned during our previous discussion. In order to 
characterize the statistical anisotropy of the CMB sky, Hajian and Souradeep [s| define 
the bipolar spherical coefficient to be, up to a sign, the expectation value of B, equation 
(JH2J)- In that respect K(l±, / 2 ; /, m), equation (JTBJ), and the bipolar spherical coefficients are 
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similar: both vanish if the CMB sky is statistically isotropic and viceversa. To reduce the 
errors of their estimates, Hajian and Souradeep also define the bipolar spectrum, which is 
quadratic in the bipolar spectrum coefficients, and hence quartic in the az m 's. We have not 
followed this path here, but have defined instead suitable estimators that are quadratic in 
the temperature fluctuations. The main difference between out approach and the one of [8( 
though is the physical focus. Whereas the authors of are concerned with the statistical 
properties of the CMB by itself, we are ultimately interested in the statistical properties of 
the primordial perturbations. Hence, our statistics K and /C (and as we shall see also A and 
A) mirror the properties of the primordial perturbations and owe its definition to a possible 
origin of statistical anisotropies in the CMB. 



C. Large Angular Scales 

Our analysis so far has been fairly general. We have not made any strong assumptions 
about the nature of the primordial perturbations, the content of the universe or the angular 
scales involved. This is why our handle on statistical anisotropies is not as strong as it 
could be. At this point we shall trade generality for constraining power. On large angular 
scales, the integral K can be explicitly evaluated, equation (|22j) . Therefore, we can use large 
angular scales to explicitly determine the power in the primordial anisotropy multipoles. 

For simplicity, let us assume that the power V\ m is scale invariant in the scales of interest, 
Vi m ~ Ai m . In a more sophisticated analysis we would try to fit for the value of the spectral 
index n. Setting accordingly n — 1 in equation (|22|) we find that the integral is proportional 
to 

T(] n= *. r[2]r[(z 2 + Z!)/2] 

1 15 2) 72 T[(h -l 2 + 3)/2]r[(/ 2 -h + 3)/2]r[(/ 2 + h + 4)/2] ' 1 ] 

which does not depend on I or m. Therefore, substituting its value into equation (|28j) we 
arrive at the following estimator for the amplitude Alm 

_ i L 1 a*i m ai+L,m+M 

LM ~ 21 + 11(1,1 + L) ^ D(l,m; L,M\l + L,m + M) ' 1 ' 

where L is again the multipole above which the power in the multipoles is cut-off. Note that 
we have actually found a set of estimators (one for each I) of the very same amplitude Alm- 
Hence, we can use this degeneracy to improve the errors of our measurements. Define the 
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averaged estimator 



AlM I -I + 



-i Imax 
t — Im.in, 



In order to improve our measurement of Alm, we would like the difference l max — l m i n to be 
as large as possible. The values of l min and l max are essentially determined by the regime 
where our large-scale approximation is valid. Therefore, following our previous discussion, 
we will set l m in = 8 and l max + L = 30 in future evaluations. Assuming again that pri- 
mordial perturbations are Gaussian and statistically isotropic we find that the mean square 
fluctuations of our estimator are 



''-max 



' \ L max l min i ^) ,_. x ' 

l — trnin 

where the variances of the individual estimators essentially follow from equation (|29|) 

/|i(0 |2\ = 1 ClCl + L V - (39) 

V LMl I I{l,l + Ly{2l + iy^D{l,m;L,M\l + Lm + My 1 ' 

So far we have concentrated on the amplitude of the primordial multipoles right below 

an eventually existing cut-off above I = L. As opposed to our previous discussion of the 

integrals K, in the case at hand it is possible to easily determine Ai m for all values of I < L 

by a recursive procedure. Imagine we know A~ lm with sufficient accuracy for I > l + l. Then, 

using equation (fTHj) we can solve for Ai m , 

(a* imi a h+lm+mi ) - Ya=i+i D (h, m 1 'J,m\h + I, m 1 + m)I(l l ,l 1 + I) A lm 
Alm ~ n77 i u — 71 ; \ • l 4U J 

Because the left hand side of the equation does not depend on m 1 , while the right hand side 
does, we could use it to construct an unbiased estimator for Ai m and compute its variance, 
just as we did before. However, since we shall not need these estimators here, we shall not 
pursue this venue any further. 

For completeness, let us mention that in full analogy with the definition of A l LM , one 
can also derive an estimator A of the amplitude Alm that transforms under an irreducible 
representation of the rotation group, 



A 



B(l,l + L;L,M) 



LM 1(1, 1 + L) Y, m i.~ l ) mD i.U m \ L, M\l + L,m + M)C(l, -m;L + L,M + M\L, M) ' 

(41) 

Its variance trivially follows from equation (JB1J); as in the case of K the variance of A is 
bigger than the one of A. As mentioned earlier, we shall not further consider the estimator 
A here. 
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FIG. 1: Plots of the real and imaginary part of K(l,l + 2; 2, 2) (data points) and its root mean 
square fluctuation in a statistically isotropic universe (continuous line). 

D. A preferred direction in the sky? 

Let us finally apply the formalism we have developed to the anisotropies we observe 
in the CMB sky. Imagine that, for whatever reason, primordial perturbations display a 
preferred direction. More precisely, suppose that primordial perturbations contain a dipolar 
component $i m . This dipole gives rise to a dipole Vim an d a quadrupole P 2 Af m the 
primordial spectrum. In the following we shall test for the existence of the quadrupole 
component. A discussion of the impact of a preferred direction on the CMB sky can be also 
found in 

Assume hence that perturbations in our universe are statistically anisotropic, and Vi m 
only vanishes for I > 3. We can then use the estimator defined in equation (fl3j) to measure 
K(l, 1+2; 2, M). If there is no power in the anisotropy multipole L = 2, then K(l, I + L; 2, M) 
vanishes. Therefore, the statistic K probes the existence of preferred direction in the pri- 
mordial perturbations. 



We have numerical 
scribed in references 



vevaluated K{1 } I + 2; 2, M) using the foreground cleaned maps de- 
27] (henceforth TOH) and [3] (henceforth LILC). Any conclusion 
derived from these maps should be interpreted with care, as it has been argued that they 
are contaminated by residual foregrounds and their noise properties are ill-understood [28]. 
However, because sky cuts in CMB maps introduce artificial anisotropies 0, at this point 
we are forced to proceed with these all-sky maps, in the hope that they accurately capture 
the real temperature fluctuations 
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FIG. 2: Plots of the real and imaginary part of (data points) and its root mean square 
fluctuation in a statistically isotropic universe (continuous line). 



For illustration, in Figure ^ we plot the real and imaginary parts of K(l, I + 2; 2, M) 
for values of I ranging between I = 2 and I = 62 and M = 2. Our choice of the upper I 
limit is dictated by our neglect of instrumental noise in the error budget. Indeed, the errors 
in WMAP's measurement of the power spectrum jsj are cosmic variance dominated (they 
contribute more than 90% of the total error) only up to a multipole of about I = 64. The 
error contours delimit the mean square fluctuations of the estimator under the assumption 
of isotropy, equation (j2HJ). In order to evaluate the latter quantities, we use the ACDM 
model that best fits the WMAP data only 

It is evident from the figure that both TOH and LILC maps lead to similar but not 
quite identical results. We can determine whether data is well-fit by statistical isotropic 
perturbations by computing the value of x 2 f° r the statistic K. In order to determine how 
X 2 is distributed, we have generated 64- 10 4 random skies drawn from a statistically isotropic 
distribution with WMAP's best fit angular power spectrum. The fraction of realizations 
P r dm{x 2 < Xo) that have an overall x 2 smaller than the one of the CMB sky in galactic 
coordinates is 97% for the TOH map and 79% for the LILC map. However, these results 
are misleading, as it turns out that the values of x 2 strongly depend on the orientation of 
the map. This is due to the non scalar nature of 

, Si ' (\K(l,l + 2;2,M)\*)' 

which is a consequence of the peculiar transformation properties of our estimator under 
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FIG. 3: The overall x 2 as a function of rotation angle. The CMB map is rotated along the x-axis 
(left) and y-axis (right) of the galactic coordinate system. Note that the plots have a period of 
180°. Chi square is invariant under rotations around the z-axis (not shown). 

rotations. In Figure E2 we plot how Xo changes as a function of the map orientation. As seen 
in the figure, the value of Xo widely varies with the orientation, so it is unclear at this point 
whether the map is consistent with statistical isotropy. 

If the CMB sky is statistically isotropic, as we rotate the temperature anisotropy map 
into all possible orientations, we expect the corresponding values of xl t° sample the dis- 
tribution of random skies -P r dm(Xo) we have considered earlier. In Figure 0] we compare the 
probability distribution function of x 2 for an isotropic sky with the distribution -P map (Xo) we 
obtain by randomly orienting the actual CMB map. At first sight, it seems that both distri- 
butions are quite similar. In order to quantify the significance of their agreement or eventual 
disagreement, let us consider the probability P that the \ 2 °f a statistically isotropic sky is 
less than the one of the randomly oriented actual CMB map, 

P = J dxl P rdm (X 2 < Xo) • ^map(Xo), (43) 

where, again, -P r dm(x 2 < Xo) i s the fraction of statistically isotropic random skies with a chi 
square less than Xo, an d -P m a P (Xo) <^Xo * s the probability that the randomly oriented actual 
CMB mas has chi square equal to xl- 

If the chi square values of a map happened to be invariant under rotations and equaled 
Xq, we would have P map (x 2 ) = Kx 2 ~ Xo), so P would simply be P rdm (x 2 < Xo)- 0n th e 
other hand, if both P r d m and P map represented the same distribution we would get P = 50%. 



19 




500 



FIG. 4: Comparison of the distribution of x 2 values of K for statistically isotropic random skies 
and random orientations of the actual CMB map. The distributions have been binned to reduce 
noise. 

Integrating equation (}4"H|) numerically we find 

P = 64% (TOH) and P = 33% (LILC), (44) 

that is, there is no significant evidence for statistical anisotropy. 

We can try to strengthen these results by considering large angular scales. In Figure [21 we 
plot the values of A% M for M = 2 and l m in = 8 < I < l max = 28. Again, as shown in Figure 
El the values of x 2 depend on the map orientation. In order to quantify up to what extent 
the maps are statistically isotropic on large angular scales, we have computed the probability 
P in equation (jIHJ), applied in this case to the chi square of A l LM . The probabilities are 

P = 51% (TOH) and P = 35% (LILC), (45) 

which agree quite well with the previous values. This is to be expected to some extent, since 
K and A are proportional to each other. However, note that the x 2 °f K an d A involve 
sums over different ranges of I. 

Finally, to confirm our conclusions, let us evaluate the estimator for the amplitude of 
an eventual quadrupole anisotropy component, equation (pT7j). In Figure El we plot how the 
corresponding x 2 varies as a function of the map orientation and compare, as in the previous 
cases, randomly simulated skies with random orientations of the actual CMB sky. In this 
case, signs of statistical anisotropics are still absent. Indeed, the probability (pl3*j) for the 
statistic A LM is 

P = 48% (TOH) and P = 43% (LILC). (46) 
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FIG. 5: In the left panel we show the overall x 2 as a function of rotation angle. The CMB map 
is rotated along the y-axis in galactic coordinates. In the right panel we plot the distribution of 
values of \ 2 for randomly generated skies and randomly chosen orientations of the LILC map. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have proposed a model-independent parametrization of the primor- 
dial spectrum that accommodates statistical anisotropies. We have determined how these 
anisotropies impact the temperature fluctuations in the cosmic microwave background, and 
what type of footprints are left behind by the former. Remarkably, the angular power spec- 
trum does not differentiate between statistically isotropic and anisotropic perturbations. As 
far as the C/'s are concerned, the angular power spectrum is as well fit by an isotropic 
primordial spectrum than by an anisotropic one. 

In order to find out whether primordial perturbations are statistically isotropic, one has 
to study the non-diagonal components of the two-point correlation function. In this paper, 
we have considered a set of statistics that measure the amplitude of the different multipoles 
in the primordial perturbations. The first set, K and A, does not transform under an 
irreducible representation of the rotation group, but has a lower variance than the second 
set, /C and A, which does possess well-defined transformation properties under rotations. In 
this paper, we have concentrated on the first set, leaving the second for future work. When 
applied to the full-sky maps of the CMB temperature fluctuations currently available, the 
former statistics do not lead to evidence for statistical anisotropy, although the analysis is 
obscured by the strong variations of our estimators with the map orientation. 

The statistical isotropy of the primordial perturbations is an observational issue. As with 
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Gaussianity, scale invariance or adiabaticity, it should be constrained rather than postulated 
or automatically implied. In this paper we have taken the first steps in that direction. If 
at the end of the journey it turns out that observational data is inconsistent with isotropic 
primordial perturbations, our current (inflationary) models for the origin of structure will 
have to be discarded and replaced by radically different ones. 
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APPENDIX A: USEFUL FORMULAE 

In this appendix we list a series of useful formulae, mainly concerning spherical harmonics. 
We follow the conventions in 31 1 (which agree with the ones of 

• Expansion of plane wave in Legendre polynomials 

exp(ik-h) = ^2(2l + l)i l ji(k)Pi(k-n). (Al) 
i 

• Addition theorem 

4-7T 

Pi(k ■ n) = ^ E Y Uk)YLifi). (A2) 

m 

• Orthonormality of spherical harmonics 

d 2 kY* m {k)Y Vm ,(k)=5 lv 5 mm ,. (A3) 

Behavior under complex conjugation 

YL = (-l) m Yi~ m . (A4) 



http : / /healpix . jpl . nasa . gov/. 
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Addition of angular momenta 

,m 1 ;l 2 ,m 2 \l,m)Yi m {k), (A5) 

where 



47T , 

Im 



{21 1 + 1) (2/2 + 1) 

D(Zi,mi;Z 2 ,ra 2 |Z,m) = x \ — 1 ^ - (h, 0; Z 2 , 0|Z, 0)(h, m X ] Z 2 , m 2 \l, m). (A6) 

The (li, mi, l 2 , rn 2 \l, m) are (real) Clebsch-Gordan coefficients. 
Symmetry properties of Clebsch-Gordan coefficients 



21 + 1 

(Zi,mi;Z 2 ,m 2 |Z,m) = (-l)' 1 '" 11 -^/ 2 ^ - 1 -mi; /, m\l 2 , ro 2 ). (A7) 

Integral of three spherical harmonics 

d 2 ^r4(A:)r /imi (fc)^ 2m2 (A;) = -£=D{h, m l5 l 2 , m 2 \l, m). (A8) 

V47T 



Note that Y^Ytm is invariant under spatial rotations. Hence, it follows from 
equations (jA8|) and (|A3|) that 



,m 2 \l,m) = (21 + 1) 5i 2 o5 m20 . (A9) 



Relation between Clebsch-Gordan coefficients and Wigner's 3-j symbol 

h h I 



(h, m i; Z 2 , mal/i, Z 2 ; Z, m> = (-l) n -' 2+m V2TTT | | . (A10) 

mi m 2 — m 
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